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1. Introduction 

In this paper we study isotropisation of Bianchi class A models with curvature when a 
minimally coupled and massive scalar field is considered. 

Locally, General Relativity (GR) with a perfect fiuid seems a good description of 
our Universe. At cosmological scale, accelerated dynamical expansion is observedjH E] 
and additional fields are required to explain it. Among them, a scalar field seems a 
good alternative although it is not the only one: higher order theories|nill] or dissipative 
fluidfnilEl also give birth to infiationary behaviour. Scalar fields are required by standard 
model for elementary particles as well as by unification theories for which, for instance, 
compactification schemes[3 [H| are considered. These last theories also give a natural 
order of magnitude for the cosmological constant [HI at early times which may be 55 to 
120 orders of magnitude bigger than its present observed value: this is the so-called 
cosmological constant problem. A solution is to consider that this "constant" varies 
across the Universe history. A massive scalar field is then an interesting possibility 
to simulate such a mechanism. All these elements show the interest of a minimally 
coupled scalar-tensor theory with a Brans-Dicke coupling function u and a potential U 
depending on the scalar field 0. 

What about the geometrical framework of this paper? FLRW models geometrically 
describe the observed homogeneity and isotropy of our Universe. However they are 
very special ones among the set of all possible models and do not allow to explain 
the observed large-scale structures. Moreover, at early times, before the decoupling 
between matter and radiation, we have no indication about Universe's geometry. Was 
it as so symmetric as the FLRW models imply? Thus, it seems interesting to generalise 
them by only keeping their spatial homogeneity property. Bianchi models describe 
anisotropic cosmological models and may allow to understand the process leading to an 
isotropic Universe. The most studied Bianchi models are those containing the FLRW 
solutions [Tfl]. i.e. the types /, V , VIIq^u and IX. We have examined the Bianchi type / 
model in [TI]. Here, we will be interested in the Bianchi class A models with curvature. 

Our goal is to look for necessary conditions allowing the isotropisation of Bianchi 
class A models with curvature when a minimally coupled and massive scalar tensor 
theory is considered. We will then deduce the common asymptotical behaviour of the 
metric functions when isotropisation is reached and compare our results with those 
obtained for the Bianchi type / model [H]. From a technical point of view, we will use the 
methods of [TT]: we will get the field equations from ADM Hamiltonian formalism |T2l fT3] 
and rewrite them with a new set of variables. Then we will look for equilibrium points 
corresponding to isotropic stable states. 

A large amount of work has been done on equilibrium states of Bianchi models. 
Wainwright, Ellis and collaborators have studied equilibrium points of homogeneous 
models for General Relativity with perfect fiuid, tilted or not and found asymptotically 
isotropic solutions. A good summary of their work is |T4]. They use Hubble-normalized 
variables to study the dynamics of Einstein field equations. The normalisation factor is 
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the Hubble parameter and the equation allowing to show that variables are normalized 
is the generalized Friedman equation. In this paper we will also consider normalized 
variables but we will use the Hamiltonian as normalisation factor. The expression for 
Hamiltonian will be the constraint from which we deduce that variables are normalized. 
More recently Barrow and KodamaflHl EE] have examined the influence of topology 
on isotropy and flatness of the Universe. They have shown that "the topology of the 
Universe can impose signiflcant restrictions upon the type of anisotropics it can sustain". 
We will not consider topology in this work but these results are really interesting from 
the point of view of relations between dynamics and topology which has also been 
examined by Ashtekar and Samuel f^. 

The plane of the paper is the following. The second section will be parsed into three 
subsections. The flrst one will be devoted to the Bianchi type // model, the second one 
to the Bianchi types VIq and VIIq and the third one to the Bianchi type VIII and IX 
models. Each of these subsections will be divided into two subsections devoted to the 
field equations and the study of the equilibrium points. We will discuss the physical 
meaning of our results in section |31 



2. Mathematical study of isotropisation for class A Bianchi models 

We begin calculating the Hamiltonian field equations. The Lagrangian of the minimally 
coupled scalar-tensor theory is given by: 

S = (167r)'i f\R- (3/2 + _ u^^^^ ^^^^^ (1) 



Although it may be more natural to redefine so that the kinetic term takes a standard 
form we prefer considering an unspecified Brans-Dicke coupling function such 

that our results be valid for any form of uj{(j)) even when it is analytically impossible to 
get 0(u;). The general form of the metric for Bianchi models is written: 

ds^ = -{N^ - NiN')dn^ + 2NidVLj + RlgijJJ (2) 

The uji are the 1-forms defining each Bianchi model. N and A^* are respectively the 
lapse and shift functions. To find the ADM Hamiltonian corresponding to the action 
(11} , we proceed as in |T2| and We rewrite the action as follows: 

S = (167r)-i j (ff''^ + n-^^ - NC^ - N,C')d^x (3) 

The Uij and H,^ are respectively the conjugate momenta of the metric functions and 
scalar field. The lapse and shift functions now play the role of Lagrange multipliers. By 
varying Q with respect to A^ and Aj, we get the constraints C° = and C* = with: 



C = ng (5) 
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We rewrite the metric functions Qij as e 2^2+2/3,^ means that fl stands for the isotropic 
part of the metric and for the anisotropic parts. Then, using Misner parameterisation 

pi = 2nUl - 2/3n6iU\ (6) 

6pij = diag{p+ + V3p^,p+ - VSp-, -2p+) (7) 

Pij = diag{P+ + V3p_,P+ - V3/?_, -2/5+) (8) 

and rewriting the Hamiltonian as H = 27rn^, from the expression ^ and the constraint 
C° = 0, we get: 

^2 J,2 



H''=pI+ pi + 12^^ + 247r'R',e-'''U + V{Q, (3.) (9) 

with p^ = nUff). The form of V{Q, P+, P-) specifies each Bianchi model and is given in 
table 1. From we derive the Hamiltonian equations: 



dH 12</)2 



dp^ (3 + 2iu)H 



(11) 



. _ dH _ dV{n,f3+,p.)/dp+ 
^+ - - 2H 

dH dV{n,P+,(3^)/d/3. 
= = 2H ^''^ 

- " 90 - ~^^3 + 2c.)i/ + ^\?> + 2u)m - (^^^ 
^_dH_dH_ ,,,206^"^"^ , 

We set A^j = and calculate that, whatever the Bianchi model, the lapse function is 
given by: 

N= ^ (16) 

Then, the relation between the time fl and the proper time t is dt = —Ndfl. 
Before starting the analysis of each Bianchi model, let us talk about some necessary 
conditions for isotropisation. By definition, the Universe isotropises if each metric 
function tends toward a common form, let us say R^. From Misner parameterisation, 
we deduce that it implies R and /3± — 0. A convenient measure of anisotropy 

is given by the a mesure of the root mean square anisotropy dp^/dt"^ + d(3_/dt^ = 
{dj3j^ / dVL^ +d(3_ / dVL^){dVL / dtY which have to decay such that isotropy occurs. Assuming 
that isotropisation is an asymptotic phenomenon arising when proper time diverges, it 
means that asymptotically dj3±/dt ^0. It is this last necessary condition only that we 
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will use in this work. It is not sufficient for isotropisation since it does not prevent (3± to 
diverge or to tend toward a big constant but it is necessary if we want that the Hubble 
factors Hi = ^g^^ be asymptotically the same, i.e. tend toward the same value dVt/dt 
in accordance with the fact that the Hubble constant is the same in any direction. Since 
the equations (fTIll) and the expression for N lead to: 

dt UnR^ ^ ' 

a stable isotropic state needs p±e^^ 0. We now look for the conditions allowing this 
limit. 

Let us assume that isotropy leads to a static Universe, i.e. tends toward a constant. 
Then p± must tend toward zero such that p±e^^ vanishes. However, from (fTB| and 
(II 2111, SI) ■ it comes that dp±/dt oc dV/d(3±e^^. Hence for the Bianchi /, VIq and VIII 
models, when /?+ and Vt tend toward some constants, the conjugate momentum diverges 
with the proper time t since 'p± tends toward a non vanishing constant and p+e^^ oo. 
Thus from the reasonable assumption that isotropy happens when t diverges, we deduce 
that isotropisation can not lead to a static Universe for these three models (it would be 
deeply unnatural that a static Universe ends for a finite value of t). For the VIIq and 
IX models, the demonstration is not so simple since then, when (3± ^ Q and VL — > const, 
dV/dp± and p± vanishes. We will show below that for these models also, isotropy 
is not possible when Q tends toward a constant. If now we assume that Q diverges, 
nothing at this stage prevents the asymptotical vanishing of p±e^^. Moreover, since /3± 
tend toward some constants for a diverging value of Q, we deduce that d[3±/dVL 
otherwise /3± would diverge with VL. Consequently, isotropisation should arise when 
VL ±oo, d(3±/dfl and pe^^ 0. These conditions are independent each other and 
of the considered Bianchi class A models. We will have to check if each of them is 
respected for each presumed isotropic equilibrium state. 

Let us compare them with isotropisation conditions defined by Collins and Hawking 
]W\ . First, as in this last paper, we have assumed that isotropisation arises when 
t — s> oo. Second, in the next sections, for each Bianchi models, we will show that 
isotropisation needs Q ^ —oo. This is the first condition that defines isotropisation 
in Collins and Hawking's paper and which implies that Universe expands indefinitely. 
Third, the fact that d(3±/dt and dp±/dQ tend toward zero satisfies their third condition 
meaning that "the anisotropy in the locally measured Hubble constant tends to zero". 
Thus from the necessary but not sufficient conditions stating that isotropisation needs 
df3±/dt — > and considering the field equations, we recover two of the four conditions 
of the Collins and Hawking definition for isotropy. It shows the physical meaning of the 
limit p±e^^ regarding isotropisation. Moreover, assuming that (3± tend toward some 
constants means that the matrix (3 whose components are the becomes a constant 
/3o which is the fourth condition defining isotropy in [T(1] . 

In what follows, we study isotropic equilibrium states of each Bianchi class A model 
with curvature. 



Isotropisation of Bianchi class A models with curvature. 



6 



2.1. The Bianchi type II model 

2.1.1. Field equations To study the equilibrium points corresponding to asymptotic 

isotropic states, we will use the following variables, specific to the Bianchi type II 
model: 

x± = p±H-^ (18) 

y = nRlVUe-^^H~^ (19) 

z = p^(f){3 + 2uj)-^/^H-^ (20) 

W = 7ri?2e-2f^+2(/3++V3/3-)^-l ^21) 

They are independent since x±, y, z and w are respectively functions of the independent 
quantities p±, 0, p^ and P±. Then, the Hamiltonian ^ yields: 

xl + xl+ 2%2 + 12^2 ^ ^2^2 ^ ^ ^22) 

We will consider this last expression as a constraint. It shows that the 5 variables 
{x±, y, z, w) are normalised. They allow us to rewrite the field equations as a first order 
equations system in the following way: 

= 72y'^x+ + 2Wx+ - 2Aw'^ (23) 
i;_ = 72y^x^ + 24w'^x^ - 2AV3w^ (24) 
y = y{Uz + 72y2 - 3 + 24^^) (25) 
z = y^{72z - 12i) + 2Aw^z (26) 
w = 2w{x+ + V3x_ + 12^2 ^ 3g^2 _ ^27) 
with i = (f)U^U^^{3 + 2uj)^^^'^. Since we want to keep lu and U undetermined, we will not 



explicit the form of i{4>). Then, the above system could not seem autonomous because 
£ = and thus it would be meaningless to look for its equilibrium point. However, it 
exists two possibilities such that it becomes autonomous. The first one is to consider i as 
a function of {x, y, z, w) rather than 0. Such considerations are often used when one look 
for exact solutions of field equations: for instance, instead of considering the potential 
as a function of the scalar field, it can be easier to find exact solutions by assuming it 
depends on the metric functions]^. However, in the general case, i may not be written 
in this way in the whole range of (p. The second possibility, which applies whatever 
i, is to consider an additional first order equation for which is derived from (fTI|l : 
(f){z,(f)). Then, the system is autonomous. However, the constraint equation shows that 
scalar field equilibrium is not necessary for isotropisation contrary to other variables: 
and can diverge whereas isotropy asymptotically occurs. Hence, we conclude that 
whatever the way we use such that the system be autonomous, equilibrium states are 
only determined by zeros of the system ()23II27|) . The same reasoning may be applied 
for each Bianchi model. In what follows, we will use (j23ll27jl to get the expressions for 
equilibrium points as some functions of £(0) and then 0(z, w) to get the asymptotical 
behaviour of the scalar field. 
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2.1.2. Isotropic equilibrium states Equilibrium points 

Several equilibrium points exist and we have to select those representing an isotropic 
equilibrium state. Immediately, we observe that the necessary conditions for isotropi- 
sation dp±/dfl 0, imply that x± —>■ near equilibrium^. Since the equations (fTIIll 
and the definitions of the variables x± will be identical for all the Bianchi models, these 
limits will be the same for each of them. Thus, we will systematically discard any equi- 
librium point with x± 7^ or which is not defined by real values. 

All the equilibrium points of the Bianchi type // model are summarised in the appendix. 
The first one is defined by {y,w) = (0,0). Let us show that it is not consistent with 
isotropy. From equations ()12II13|) . we derive that p+ — v^p_ = Po,Po being an integration 
constant. Then, using (fTTHl . it comes that /3+ — -\/3/3_ = poH~^. Thus, isotropisation 
needs a diverging Hamiltonian but for the zero measure case with po = 0. If now we 
consider (f27|l near {y,w) = (0,0), we deduce that w behaves as e~^^ and vanishes when 
Q +00. Introducing this expression in (pT| . we derive that H should asymptotically 
behave as e^('^++^^-) and thus should tend toward a constant value when Universe ap- 
proaches isotropy. This contradict the fact that the Hamiltonian have to diverge and 
thus the set of points {y,w) = (0,0) is not compatible with isotropisation. 
The second set of points is not real and such that x± 7^ 0. Hence, it does not 
correspond to an isotropic equilibrium state and we discard it. The third set of 
points is such that x± = and respects the constraint. It writes {x+,x_,y, z,w) = 
(0, 0, ±V3 — P{Qy/2)~^, i/Q, 0) and is real if £^ < 3, implying that isotropisation is not 
possible if this last quantity diverges§. Indeed since asymptotically z behaves as £/6, it 
means that an isotropic stable state needs i to tend to a constant value with no limit 
cycle otherwise z -/^ when Q —00. This is the only set of points representing 
an isotropic stable state and the only one we will consider below. The fourth set of 
points does not represent a stable isotropic state except if £ — > 1. Then it tends toward 
the previous one. However, we will see below that the value i = 1 does not allow for 
isotropisation. Thus we discard it. 

Monotonic functions 

We rewrite the equation (fT5|l with the normalised variables: 

H = -H{72y'^ + 24w^) (28) 

Hence the Hamiltonian is a monotonic function of fl with a constant sign. Then, from 
the lapse function expression (fT6|l . we deduce that is a monotonic function of the 

I Note that we have then, near equilibrium, x± ^ and x± 0. This is a consequence of the field 
equations and values of the equilibrium points near isotropy. It means that x± should be integrable in 
the Lebesgue sense in the neighbourhood of equilibrium. We will see that it is actually the case when 
we will calculate the asymptotical behaviours of x± . 

§ We will not take into account the case for which would have a chaotic behaviour such that it stays 
smaller than 3. 
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proper time t. Therefore, if the Hamiltonian is initially positive (negative), f2 —oo 
corresponds to late time(early time). We will not consider the case Q — > +00 since we 
will show below that it does not lead to isotropy. We conclude that late times isotropisa- 
tion initially needs H > 0: this is the only necessary initial condition for this behaviour. 
Moreover, the Hamiltonian being of constant sign, it is the same for the variables y and 
w. 

Asymptotic behaviours 

We wish to evaluate the behaviours of some quantities in the neighbourhood of the 
equilibrium points {x+, X-,y, z,w) = (0, 0, ±V3 — £^(6\/2)~^, i/6, 0), i.e. when we ap- 
proach isotropy in — > —00. Approximating (f27ll near equilibrium by w{l — i^), we 
find that asymptotically w behaves as e^^"^^)^. 

From this last expression and approximating (f23|l by (3 — £^)x+ — 24w^, we deduce 
that x+ behaves as the sum of two terms e^^^"^^)^ and e*^^"^^)^. Since isotropy needs 
x+ — > and £^ < 3, we derive it only occurs if Q —00 and £^ < 1. These two 
limits are in accordance with the vanishing of x± and w which is necessary to reach the 
equilibrium isotropic state. Consequently x± asymptotically behave as e^*^^~^^)^. Let 
us note that the two limits x± — and Q — > —00, necessary for isotropisation, are 
compatible and justify our assumption that it takes place for a diverging value of t: 
when — > —00 and the Universe isotropises, it is expanding and there is no physical 
meaning to consider that it ends for a finite value of the proper time. 
Let us show that the value £^ = 1 has to be discarded. If 1 — £^ — faster than w 
tends toward a non vanishing constant and is not compatible with isotropy. If 1 — £^ — > 
slower than from (fT5|) we deduce that near equilibrium H e~^^. Then from (fT2|l . 
it comes that p+ — > e~^^ and it follows from (fT8|) that x+ tends toward a non vanishing 
constant. Hence, the limit £^ — > 1 is not compatible with isotropy. The above reasoning 
concerning the case for which 1 — £^ — > faster than fl~^ will stay valid for any Bianchi 
model. However, when 1 — i"^ —>■ slower than it will be valid only for Bianchi 

type //, VIq and VIII models. For Bianchi type VIIq and IX models, the situation is 
different because when P± —>■ during isotropisation, p± varies slower than e~^^ since 
dV/dp±e^^ 0. However, we will show below that near equilibrium, the value i = 1 
can also be excluded. 

Now it is possible to show that p±e^^ vanishes when Q ^ —00. Writing p±/II as 
a function of x± and w and using their asymptotical behaviours, we calculate that 
p±/p± tends toward the constant —(1 + £^). Consequently p±e^^ ^{2-e^)n _^ g ^j^gj^ 
Q —00 and the necessary but not sufficient condition for isotropisation is respected. 
Moreover, from (f28|l . we calculate that asymptotically IIII~^ tends toward £^ — 3. Thus 
H —>■ e^^^"^)^ and diverges since < 1: therefore, although determined independently, 
the asymptotical behaviours of p± and H agreed with these of x± = p±II~^. 

Concerning the scalar field, we can find a differential equation whose solution 
asymptotically behaves in the same way as (f) when — — 00 by expressing equation 
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with the normalised variables and considering its asymptotical limit near equilibrium. 
It comes: 

) ^ Y 29 
^ (3 + 2cj)f/ ^ ' 

This important equation allows us to get the asymptotical behaviour of (p near 
equilibrium and consequently that of i. 

From the asymptotical behaviour of H and the expression (|TH|l for the lapse 
function, it is possible to get the isotropic part of the metric, e~^, as a function of 
the proper time. If i'^ tends toward a non vanishing constant, then — ^ ^ If 
tends to faster than (— fi)^^, behaves like an exponential. Let us show that it 
is always the case. Equation (f29|l can be rewritten as d(j)/dQ = 2l'^U{U^)~'^ and then 
U (X e'^-f ^ '^^ when f2 — — oo. Introducing this expression and the expression for H into 
(fT9|) . we get y oc ^-^'^^+1 ^'^'^^ _ Thus, if i"^ vanishes slower than (— y diverges or 
tends to instead of a non vanishing constant and there is no equilibrium. 

Hence, when an isotropic equilibrium state is reached with £ — > 0, always van- 
ishes faster than (— and Universe always tends toward a De Sitter one. This proof 
for £ = relies on the asymptotic form of 0, H and the definition of y. It will be valid 
for all the Bianchi models since we will see that these 3 quantities keep the same forms 
each time the same set of equilibrium points is considered. 

From (fT9|l and the asymptotical forms for fl{t), we deduce that the potential behaves 
as when i tends toward a non vanishing constant, or as a non vanishing constant 
otherwise (the same behaviour held when we consider the Bianchi type / modelfTI]). 
Concerning the 3-curvature which can be expressed as R^^^ = e2^+^(^++^^-), it is ob- 
vious that near isotropy, it tends to zero, showing that the Universe becomes spatially 
fiat. 

Let us note the importance of the potential for isotropisation of the Bianchi type // 
model. If we consider U = 0, we get from the field equations (fT2|) and (fT5|l that 
H = + Po and thus /3+ = 1 —poH^^, po being an integration constant. It follows that 
does not asymptotically vanish and isotropy can not be reached. 

Partial equilibrium\\ 

Above we have defined the isotropic equilibrium states such that all the variables reach 
equilibrium. However, this last statement is not mandatory: all of them have not to 
reach equilibrium such that x-t — ^ in a stable way. In this case, as — > —oo, the 
concerned variables would stay bounded but their derivatives would not asymptotically 
vanish: they should oscillate indefinitely (around a constant or not) and thus, their 
derivatives should also oscillate around zero without beeing damped. What are the 
variables able to behave in this way? 

We can reasonably assume that it is not x±, otherwise it would mean that the variation 
II I thank one of the referees criticisms from which this section is inspired. 
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of the Hubble constant would be anisotropic. Since x± and x± have to vanish, it implies 
that w also asymptotically vanishes and consequently, w — > 0. Finally the only variables 
whose equilibrium is not necessary to isotropy and whose derivatives could be oscillating 
are y and z. Under these assumptions, what about C behaviour? 

Equation (pHl) writes asymptotically z = y'^{72z — 12i). Since we have assumed that 
z was oscillating, it shows that £ can not tend to a constant, diverge monotonically or 
diverge oscillatory if the oscillations are not large enought(£ = fi^/^ + cosfi as instance) 
otherwise the sign of z would be constant. Consequently, only oscillatory £ with suffi- 
ciently large amplitudes and not tending to a constant may allow oscillations of z too 
{£ = flsinfl or £ = n cos Q with n a constant larger than the largest amplitude of 6z 
as instance). In this case, the results of the previous sections do not apply since all the 
variables do not reach equilibrium but an isotropic equilibrium state eventually occurs 
with w) only reaching equilibrium. Then, since £ can be a regular functions as well 
as having a chaotic behaviour, it is not possible to give more characteristics about it or 
the way x± would reach equilibrium. Hence, the main result of this subsection is mainly 
a limitation of the previous subsections results which will be valid for all the Bianchi 
class A models. 

Let us examine the folowing example for the Bianchi type II model: 
Since x± and w vanish, asymptotically the constraint is 2%^ + 12z'^ = 1. In this limit, 
considering that ipHl) is essentially equivalent to the constraint equation under (^^1) . (^^1) 
is the only nontrivial equation. In this asymptotic reduction, (|25|1 can be regarded as 
the equation for £ in terms of y, which can be written as: 



2 + Vy/T 



where v = 2(36?/^ — 1) and 6z = \/l — v. For any function for v{Q), the equation 

. 120 

= , z, 
V3 + 2cj 

determines as a function of fl for a given u. Then, the definition of £, 





V3 + 2w U 

determines f/ as a function of 0, provided that (p — Q relation is invertible. Here, if 
V is bounded by a positive constant from below and if 1 — f is non-negative, one can 
easily check that x± and w have required asymptotic behavior. Thus, the only possible 
constraint on v{fl) is the invertibility of the cj) — fl relation. For example, the choice: 

= t^o + ^ sin^ 

satisfies this condition, if Vq is a constant in the range < fo < 1. However, for this 
choice, £, y and i diverge oscillatorily as — — oo, although the isotropization condition 
x± ^ as Q — oo is satisfied and {y, z) are bounded. 
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2.2. The Bianchi type VIq and VIIq models 

The results are similar to those of Bianchi type // model. 

2.2.1. Field equations We will use the following variables for both Bianchi type VIq 



and VIIq models: 

x± = P±H-^ (30) 

y = TcRle-^'^U'/'H-' (31) 

z = p^(f){3 + 2n)-^/^H-^ (32) 

They are independent since they respectively depend on p±, 0, p^ and a combining of 
/?±. The definitions of a;±, y and z are the same as these of the Bianchi type // model. 
The Hamiltonian is written: 

xl + xl+ 2%2 + 12^2 ^ ^2{w+ ± w_Y = 1 (34) 
and the field equations become: 

x+ = 72y'^x+ + 2A{x+-l){w^±w+)^ (35) 

X- = 72?/\_ + 24x_ {w- ± w+f + 24V3{w'i - wl) (36) 

y = y{6iz + 72y^ -3 + 24:{w^±w+f) (37) 

z = y^{72z-12£) + 2Az{w.±w+f (38) 

w+ = 2w+ [x+ + VSx- + 12{w^ ± w+)^ + 36y^ - l] (39) 

w- = 2w- \x+ - V3x_ + 12(«;_ ± w+)^ + 36y^ - l] (40) 



The ± symbols in equations ()34ll40p correspond respectively to the Bianchi type VIq and 
VIIq models. For the first model, the constraint shows that the variables are normalised 
since w+ and W- are positive. For the second one, because of the minus sign, both w+ 
and w_ could diverge. Then the constraint will be respected only if the sum w+ — 
tends toward a constant. We will show below that isotropy is only compatible with finite 
values of w±. Consequently, the isotropic states we are looking for are reached for some 
bounded values of the variables whatever the Bianchi type VIq or VIIq models. Now, 
we can also show that isotropisation of Bianchi type VIq and VIIo models may not 
arise for a finite value of f2. Indeed, if f2 tends toward a constant when the proper time 
diverges, dfl/dt vanishes and thus, from (fT?)| . it comes that H vanishes. But then, w± 
should diverge which is not compatible with the equilibrium as shown above. Hence, for 
the Bianchi type VIq and VIIq models, isotropisation does not lead to a static Universe. 



2.2.2. Isotropic equilibrium states Equilibrium points 
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All the equilibrium points with finite values of and w_ are referenced in 
the appendix. Following the same reasoning as for Bianchi type // model, 
the only set of equilibrium points compatible with an isotropic stable state is 
{x+,x-,y,z,w+,w-) = (0, 0, ±V3 - i^{6V2)-\i/6, 0, 0), implying that < 3. It 
is equivalent to the points found for Bianchi type II model and, for the same rea- 
sons, i have to tend to a constant with no limit cycle. Another interesting set 
of points is given by {x+,X-,y, z,w±,w±) = ((^^ - + 8)"\ ±^3(^2 - + 
8)-\ ±V12 - 3& [2{e + 8)]'^ , U{2e + 16)-\ 0, + 5£2 - 4 [2{e + 8)]"^). How- 

ever, in this case x± ^ only if — ^ 1 and we recover the values of the previous set 
of points for this particular limit of £^ which does not allow isotropisation as it will be 
shown below. Other sets exist but are complex valued and can thus be discarded. At 
last, as for Bianchi type // model, let us show that the set of equilibrium points such 
that = (0,0,0) implies that x± do not vanish. 

In this case we deduce from ()39ll40p that w± behave as e"^^ and thus isotropy would 
arise when VL — > +oo. Then, using the definition for w±, we derive that H should 
be a constant near isotropy and, considering equations ()12II13|I . we find that asymptoti- 
cally p± should vary as Po±e~^^ (po± being some constants) for Bianchi type VIq or even 
slower for the VIIq type since l3± might tend toward some vanishing constants. It would 
follow that p± pi-t, pi± being some integration constants. However, in this case (3± 
would tend asymptotically toward the constants pi±H~^ and isotropisation could not 
occur for a diverging value of VL. Hence, (y, w+,i(7_) = (0,0,0) is not compatible with 
isotropisation except for the special case of zero measure pi± = 0. 



Monotonic functions 



What about monotonic functions? We can rewrite equation (|T5|) as follows: 



H 



-H 



72y^ + 24:{w+ ± w_ 



(41) 



As for the Bianchi type II model, (|4T| shows that the Hamiltonian is a monotonic 
function of constant sign which determines if isotropisation occurs at early or late times 
depending if the Hamiltonian is initially negative or positive. Moreover, it follows that 
y and w± are of constant sign. 



Asymptotic behaviours 

Making the same approximation as for subsection I2.1.2| we find that near an isotropic 
equilibrium state, w± — > e*^^"^^)^. Then, assuming that isotropisation arises for Q — > —oo 
as we will show it below, it follows from (jlT|l that H e~^^ when £ — > 1 and then, 
near equilibrium, w± should tend toward some non vanishing constants. Thus the value 
i = 1 does not agree with isotropisation. 

Concerning x±, from ()35II36|) . we deduce that they asymptotically behave as the sum 
of two exponentials e^^^"^^)^ and e^^^^^^^, showing again that these quantities will tend 
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toward zero only if 17 — oo and P < 1. These two limits allow vanishing oiw±, which 
is necessary to reach the equilibrium states. 

All these elements show that near isotropy w± is bounded. Indeed, for (x, y, z) to reach 
equilibrium, one only needs that w+ — W- — > whatever the particular asymptotical 
behaviours oiw±. Then only using this last limit, we have recovered the asymptotical 
behaviours for x± which implies that < 1. From ()39ll4np and these last expressions 
and condition, it is then possible to get the particular behaviours of w±, which show 
that, near isotropy, these variables always vanish and are then bounded. 
Again, we are able to calculate that p±e^^ tends to as e(^~^^)^, thus filling the necessary 
but not sufficient conditions for isotropisation we defined above. From equation (|1T]1 we 
derive that asymptotically H behaves as e^^^"^)^. The form of equation (fTT|) for being 
unchanged whatever the Bianchi model as well as those of I and z near equilibrium, we 
recover the same differential equation as giving the asymptotical behaviour for the 
scalar field. 

Since H and when ^ — oo, and N and y have the same forms as for the Bianchi 
type II model, we find the same behaviours for and [/ as a function of the proper 
time depending if ^ tends or not toward a vanishing constant. Again, the 3-curvature 
tends to zero since (3± become constant when VL diverges negatively. 



2. 3. The Bianchi type VIII and IX models 

2.3.1. Field equations We will use the following variables: 

x± = p±H~^ 

y = T^Rle-^'^U'/'H-' 

z = p^(f){3 + 2n)-^/'^H-^ 



The variables x±, y and z are the same as those defined for the Bianchi type // model. 
Wp and Wm are not independent because both are related to Near isotropy, we will 
have Wm oc Wp oc e~^^II~^. w_ is positive. The constraint equation is written: 

xl + x'i+ 2%2 + 12z^ + I2{wl{l + wt) ± 2w-{wmWpf^^{l + w!) + 
wl{wl-2wl)]{w'^wpr' = l 



and for the field equations it comes: 



x+ = 72y'x+ + 24{w;J(x+ - 1)(1 + wt) ±w^{l + 2x+){wmWpf/\l + w'i) 



+ wt 

x_ = 



{2 + x^)wl~2{x+-l)wl\}{wlwp) 
72y'^x^ + 24:{wl 



-1 



(42) 



w"^ (x^ 



V3) + X_ + VS) \ ± W^{w.mWpf^^[w^_ 
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(-V3 + 2x_) + (VS + 2x_)] + w^^x^wl - 2wl)}iw^_Wp)-^ (43) 
y = y{6iz + 72y^ - 3 + 24[u;J(l + wt) ± 2(u7„u7p)3/2^_(l + w^) + 

w'_{wl-2wl)]{w'_Wp)-'} (44) 



z = y\72z - 12£) + 2Az[w;{l + wt) ± 2{w^Wpy''w_{l + wi) + 
wliwl-2wl)]iwlwp)-' (45) 
Wp = Wp{-2 + 2x+ + 72y^ + 24[w|(l + wt) ± 2u;_(w^Wp)=^/2(x ^ ^2 ^ 
+ njl{wl~2wl)\{wlwp)-'} (46) 

= ^i;„^{-2 - 2x+ + 72?/2 + 24[u7j(l + wt) ± 2w_(u;™Wp)^/2(l + wt) 
+ wl{wl-2wl)]{wlwj;)-'} (47) 
w;_ = 2VIW-X- (48) 



± being related respectively to the Bianchi type VIII and IX models. For sake 
of completeness, we have written differential equations for each variable Wp and Wm- 
However, they are equivalent. The constraint shows that the variables are not necessarily 
normalised: if one of them diverges near isotropy, it have to be counterbalanced by the 
divergence of Wp and Wm- Thus, if we show that isotropy does not arise for unbounded 
values of Wp and Wm, it will mean that it only happens for some finite values of the 
variables. 

To reach this goal, we will write that Wp Wm — ^ w and W- 1. This is justified 
because isotropy needs P± ^ and we will see below that an isotropic equilibrium state 
effectively implies W- ^ 1. In this case, for Bianchi type VIII model, all the variables in 
the constraint are positives and thus bounded. Concerning the Bianchi type IX model, 
let us assume that w diverges. Then, putting x± = 0, from the constraint we have 
asymptotically 3w^ 2y'^ + z"^ — 1/12 and from the equation for w, 3w^ 3y^ — 1/12, 
implying that asymptotically z'^ y"^ and diverges as w"^. However, with these limits 
we get from the equations for y and z that y — > 6iz^ — 3z and i —12iz'^ + 2z. Then 
equilibrium for y and z can only be reached if z ^ which is in disagreement with the 
assumption on the divergence of w. Hence, an isotropic equilibrium state is not possible 
if Wp and Wm diverge. It follows, in the same way as for the Bianchi type VIq and VI Iq 
models, that isotropisation of Bianchi type VIII and IX models for a finite value of fl 
is impossible. 

We can also show that Wp and Wm may not tend toward some non vanishing constants. 
Let us assume that it is actually the case and consider 2 constants w and a such that 
Wp ^ w and Wm ctw. We introduce these limits in equations (|42ll43p with x± = 
and get respectively: 

x+ = -2Aw^{l + w_a^/^(l + wt) - 2wtil + a^) + wt)wZ'^ (49) 
x_ = -24:V3w'^{wl - 1)(1 - + wt)wZ^ (50) 
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Then, for Bianchi type VIII model, we derive that equilibrium for x± will be reached 
only if a tends toward a complex value (—1)^/^ or/and W- is negative, which is im- 
possible. For the Bianchi type IX model, equilibrium for x± may be reached if 
Wp Wm (i.e. P± 0) and w_ — > 1. Then, looking for the equilibrium points, 
the only ones which may be real and such that Wp and Wm be different from are 
{x+,x_,y,z,Wp,Wm,w_) = (0,0,±(6^)-\(6£)-\±(l - £2)i/2(Q^)-i)_ xhey check the 
constraint equation and are real if £^ < 1. Then we calculate that Wp and Wm behave 



like± 



£2)(l-e P — -+36f 



1/2 



and thus reach equilibrium when Cl +oo. 
Meantime, starting from this last expression and introducing it in the equation for x+, it 
comes that x+ tends toward a complex value when f2 +oo and thus these equilibrium 
points are excluded. Numerical simulations seem to confirm the absence of equilibrium 
for these values of {x+, X-,y, z, Wp, Wm, W-). 



2.3.2. Isotropic equilibrium states Equilibrium points 

To find the equilibrium points we have to consider the equations ()42II45|) . (|i8|) and 
one of the equations (HHll or (jTZjl since both Wm and Wp depend on j3+. However the 
solutions of the equations system thus defined can not be easily calculated. Conse- 
quently, we will only take into account the solutions such that {x±,Wp,Wm) = (0,0,0) 
and which are compatible with isotropy. Then the solutions reduce to the set of equilib- 
rium points {x+, x_,y, z,Wp,Wm,w^) = (0, 0, ±y/3 — P{6\/2)~^, i/6, 0, 0, 1) which will 
be real if £^ < 3 and respect the constraint equation. It is equivalent to the sets found 
for the previous models and once again, i have to tend to a constant with no limit cycle 
such that equilibrium be reached. Note that it is such that /3_ — > since 1. 



Monotonic functions 



We can rewrite ()15|1 in the following form: 

H = -H[72y' + 24(±2 ^^ ± 2w'J'w'jl'w. - 2wl + ^ + 

wlwl + ^)] (51) 

Wp 

(52) 

We immediately see that it is not a monotonic function and that its sign is indefinite. 
Thus Q is not a monotonic function of t and it is not possible to determine if isotropi- 
sation, corresponding to i7 ^ —oo, arises at early or late proper times. 



Asymptotic behaviours 



Near equilibrium, it is possible to approximate equation by Wp{l — ^'^) implying that 
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Wp tends toward e^^"^^)^. The same conclusion arises for Wm- In the same way as the 
previous subsection, one can show that the value £^ = 1 is not agreed with isotropisa- 
tion. Introducing asymptotical expressions for Wp and Wm in the equations (j42ll43|) . we 
find that x± behave as the sum of two exponentials, e^*^^~^^^^ and e^^"^^)^. Thus, once 
again, isotropisation needs fl — > — oo and £^ < 1 implying that x± behave as e^(^~^^)^. 
As for Bianchi type // model, it is possible to show that p±e^^ — > e*^^"^^^^ and thus van- 
ish. Moreover H behaves again as e^^^"^)^. The Hamiltonian equation (fTT| for being 
independent of considered Bianchi model and the definition and asymptotical value of 
z being the same as for Bianchi type // model, we find the same differential equation 
(f^ . giving asymptotically the behaviour of the scalar field. Anew, since H and 4> when 
Q — oo, and and y have the same forms as for Bianchi type II model, the dis- 
cussion about the forms of the metric functions near equilibrium is the same and they 
behave as power or exponential law of the proper time depending on the asymptotical 
value of i. 

3. Discussion 

We have found some necessary conditions for isotropisation of Bianchi class A models 
with curvature for a minimally coupled scalar tensor theory. We have seen that the Uni- 
verse has to expand {fl —oo), justifying the assumption that t should be diverging, 
and that the ratio between the conjugate momentum and the Hamiltonian should van- 
ish. Our results do not concern the class of theories for which i prevents the equilibrium 
of z and y. As shown in subsection I2.1.2| such £ should be oscillating with significant 
amplitude and not tending to a constant. Hence, they concern the i which tend to a 
constant, diverge monotonically or even with negligible oscillations. In these cases, our 
main result is: 

A necessary condition for isotropisation of Bianchi class A models with curvature for 
General Relativity plus a massive scalar field, whatever the considered Brans-Dicke cou- 
pling function and potential, will be that (f)U^U~^{?) + 2uj)~^^'^ tends to a constant £ such 
that < 1. For Bianchi type II , VIq and VIIq models, it arises at late times if the 
Hamiltonian is positive, at early times otherwise. For Bianchi type VIII and IX mod- 
els, the time of isotropisation is undetermined. If isotropisation arises with £ the 
metric functions tend toward a power inflationary law t^ ^ and the potential vanishes as 
t~^. If it arises as i = 0, the Universe tends toward a De Sitter model and the potential 
to a constant. In any case, isotropisation requires late time accelerated expansion and 
the Universe becomes spatially flat. 

Necessary condition for isotropisation determines an asymptotical limit that the scalar 
field have to respect. It can be compared to the limit required such that scalar tensor 
theories be compatible with solar system tests when ?7 = 0, i.e. ^ oo and uj^j^uj'^ 0. 
To evaluate £, we need to know the asymptotical behaviour of the scalar field. It comes: 
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The value of the scalar field when the Universe reaches an equilibrium isotropic state 
with an asymptotically constant I is the value of the function (j) defined by (p = 
202[/^(3 + 2uj)-^U-^ when n -oo. 

Although Bianchi type IX model contains the closed FLRW solutions, when the Uni- 
verse isotropises and we consider a minimally coupled and massive scalar field, it is 
infinitely expanding. Moreover, the common late time attractor of all the isotropising 
solutions is not oscillating. This fact may seem astonishing for the Bianchi type VIII 
and IX models. However, in it has been observed that despite the mixmaster 
behaviour of Bianchi type VIII model at early time, its late time behaviour can be 
non oscillating. If we compare the results got when no curvature is present [TT] with 
those of this paper, few differences appear. The asymptotical behaviours of the scalar 
field and isotropic part of the metric are the same in both papers, partly because the 
Hamiltonian and lapse function asymptotically behave in the same way. However the 
variations of the functions describing the anisotropy, f3±, are different since the conju- 
gate momenta are not constant in presence of curvature. The fundamental difference 
comes from the interval of i allowing for isotropy. For the Bianchi type / model, it 
was £^ < 3 and decelerated dynamics was possible. For the models with curvature, 
we have £^ < 1, implying that isotropisation requires late time accelerated expansion. 
This is due to the presence of curvature which reduces the interval of values of i related 
to the Bianchi type I model. Hence, late time accelerated expansion finds a natural 
explanation through the fact that our Universe is isotropic. Other problems are nat- 
urally solved by isotropisation: asymptotically, the 3-curvature vanishes thus solving 
the fiatness problem. It comes from the fact that during isotropisation P± tend toward 
a constant whereas Q ^ —oo. In the same way, the small value of the cosmological 
constant could be explained by the fact that when Universe isotropises and £ does not 
vanish, the potential, which mimics a dynamical cosmological constant vanishes. If £ 
vanishes, the potential tends to a non vanishing constant but it is not necessary small 
except if we fine tune it. 

To complete this study, let us write few words about Bianchi class B models. Their 
Hamiltonian formulation is different from those of class A and has been studied in [22] . 
It needs to redefine the divergence theorem in a non-coordinated basis. Then the Bianchi 
type V Hamiltonian writes as the one of Bianchi type / with an additional constraint 
p+ = 0. Consequently for Bianchi type V model which contains the solutions of open 
FLRW model, isotropisation follows the same rules as those of Bianchi type / model 
described in |TT|. The nature of the other class B Hamiltonians is totally different and 
will not be considered here. 

Let us examine some results usually considered in the literature. 

The "No Hair theorem" of Wald[23] states that General Relativity with a scalar 
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field and a cosmological constant isotropises toward a De Sitter Universe. Here, as for 
the Bianchi type / model, when £ and if the minimally coupled scalar tensor theory 
isotropises^, it will tend toward a De Sitter Universe. This generalise the "No Hair 
theorem" which takes into account only the case U = ate for which i = 0. 

It is really shocking that it exists only one set of equilibrium points shared by all 
the Bianchi models and representing the only possible isotropic stable equilibrium state. 
However, despite a careful analysis we have not found any additional points with such 
properties. One way to check if this statement is true is to select some special forms of 
U and u and then to verify if the conditions for isotropisation of the theory thus defined 
and the asymptotical value of are in agreement with our results. It can be easily 
done with the theory defined by an exponential potential U = e^'^ and a Brans-Dicke 
coupling function ^/3 + 2uj(f)~'^ = ^J2 whose isotropisation has been extensively studied 
in the literature using different methods. In this case P' = k'^/2. We have collected 
the conclusions of different papers and have compared them with ours. In 1231 [^Hl, it 
is shown that isotropisation arises at late time when fc^ < 2 (except the contracting 
Bianchi type IX models) and lead to a De Sitter Universe when = or to a power 
law of the form ^ for the metric functions otherwise. If /c^ > 2, the Bianchi type /, 
V, VII and IX models might isotropise at late times. Concerning the Bianchi type / 
model, we have shown in [TT| that a necessary condition for isotropisation will be fc^ < 6 
but it was impossible for larger values. For the models of class A with curvature, from 
the present paper we deduce that isotropisation is possible only when < 2 and always 
comes with late time accelerated expansion. The asymptotic behaviour of the metric 
functions is in accordance with that of [23112^1. A difference is that Bianchi type VIIq 
and IX should not isotropise if k'^ > 2. Concerning the Bianchi type VIq model, our 
results agree with these of [22]. For the Bianchi type V model, they are the same as 
these of the Bianchi type / model in accordance with [23]. Hence, concerning the special 
case of an exponential potential, there are few differences between our results and those 
of others papers. It seems to confirm the presence of a unique set of equilibrium points 
shared by all the Bianchi class A models and representing an isotropic equilibrium state. 
Of course, the case of an exponential potential could be a particular one and thus other 
types of potentials should be studied to check the results of the present paper. Note 
that, isotropic state is not the only possible late time equilibrium state. As written 
above or shown in the appendix, other ones exist, for instance with x-t 7^ 0, but they 
do not correspond to an isotropic Universe. 

To conclude. Universe isotropisation requires late time accelerated expansion 
because of the curvature. Then, it becomes fiat and the potential vanishes as 
or tends toward a constant. These features fit well with the current observations and 
leave the door open to geometrical and physical generalisations of standard cosmological 
framework. In a next work, we will take into account the presence of a perfect fiuid. 

% Do not forget that £ tending toward a constant is a necessary but not sufficient condition for 
isotropisation. 
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4. Appendix 

In this appendix, we present all the equilibrium points of Bianchi type //, VIq and VIIq 
models. 



Bianchi type II model : 
. (y,u;) = (0,0) 

• = (1, V3, 0,0,±V2) 

. xi, z, w) = (0, 0, ^/6, 0) 

\x+,xi,y,^,w) l£2_,_8, V'J^2_,_8, =c 2(^2_,_8) , 2(£2+8)' ^ 2(£2+8) ) 

Bianchi type VIq and VIIq models : 

. {y,w+,w^) = (0,0,0) 

• x-,y,w+,w^) = {1,0,0, w+,w+) 

• x_,y, z, w+, w_) = (1, -VS, 0, 0, 0, ±i/2) 

• x_,y, z, w+, w_) = (1, 0, 0, ±i/2, 0) 

. x_, y, z, w+, w.) = (0, 0, i/6, 0, 0) 



yx+,x-,y,z,uj+,uj-) — \^2j^^, V o^a+g , =c 2(^2_,_8) ' 2(^2+8) ' ^' ^ 2{£2+8) 

fr T 1, 7 m m \ - a/^^'-^ + v^12^ ^ v/(€2_i)(^2_4) 

x_, (/, — l,£2_,_8' V -Ji^qig, =^3{F+8)"' 2(F+8) ' ^ 2{F+8) ' 
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Table 1. Form of V{n, f3+, j3^) for each Bianchi model. 

Bianchi type V{n,l3+,I3-) 

II 127r2i?4e4(-n+/3++y3/3_) 

VIo, VIIq 247r2i?4e-4^2+4/5+ (cosh4\/3/3„ ± 1) 
VIII, IX 247r2i?4e-40[g4/3+(gQgj^4^^_ _ 

1/26-8/3+ _^ 26-2/3+ cosh2\^/3^] 



